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ABSTRACT 

The orbital evolution of about 20000 Jupiter-crossing objects and 1500 resonant asteroids under 
the gravitational influence of planets was investigated. The rate of their collisions with the 
terrestrial planets was estimated by computing the probabilities of collisions based on random- 
phase approximations and the orbital elements sampled with a 500 yr step. The Bulirsh-Stoer 
and a symplectic orbit integrators gave similar results for orbital evolution, but sometimes gave 
different collision probabilities with the Sun. For orbits close to that of Comet 2P, the mean 
collision probabilities of Jupiter-crossing objects with the terrestrial planets were greater by two 
orders of magnitude than for some other comets. For initial orbital elements close to those of 
Comets 2P, lOP, 44P and 113P, a few objects (~0.1%) got Earth-crossing orbits with semi-major 
axes a<2 AU and moved in such orbits for more than 1 Myr (up to tens or even hundreds of 
Myrs). Some of them even got inner-Earth orbits (i.e., with aphelion distance (5<0.983 AU) and 
Aten orbits. Most former trans-Neptunian objects that have typical near-Earth object orbits 
moved in such orbits for millions of years (if they did not disintegrate into mini-comets), so 
during most of this time they were extinct comets. 



INTRODUCTION 

The orbits of more than 70,000 main-belt asteroids, 1000 near-Earth objects (NEOs), 670 
trans-Neptunian objects (TNOs), and 1000 comets are known. Most of the small bodies are 
located in the main asteroid and trans-Neptunian (Edgeworth-Kuiper) belts and in the Gort 
cloud. These belts and the cloud are considered to be the main sources of the objects that 
could collide with the Earth. About 0.4% of the encounters within 0.2 AU of the Earth are 
from periodic comets (http://cfa-www.harvard.edu/iau/lists/CloseApp.htmll, and 6 out of 20 
recent approaches of comets with the Earth within 0.102 AU were due to periodic comets 
(http:/ /cfa- www. harvard.edu/iau/lists/ClosestComets. html ). So the fraction of close encounters 
with the Earth due to active comets is ~1%. Reviews of the asteroid and comet hazard were 
given in [l]-[3]. Many scientists [3]- [5] believe that asteroids are the main source of NEOs (i.e. 
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objects with perihelion distance q<1.3 AU). Bottkc ct ah [3] considered that there are 200ibl40 
km-sized Jupiter-family comets at q<l.3 AU, with ~80% of them being extinct comets. 

Duncan et al. [6] and Kuchner [7] investigated the migration of TNOs to Neptune's orbit, and 
Levison and Duncan [8] studied the migration from Neptune's orbit to Jupiter's orbit. Ipatov 
and Hahn [9] considered the migration of 48 Jupiter-crossing objects (JCOs) with initial orbits 
close to the orbit of Comet P/1996 R2 and found that on average such objects spend 5000 yr 
in orbits which cross both the orbits of Jupiter and Earth. Using these results and additional 
orbit integrations, and assuming that there are 5 x lO'' 1-km TNOs with 30<a<50 AU [10], 
Ipatov [1],[2],[11] found that about 10^ 1-km former TNOs are Jupitcr-crossers now and 10-20% 
or more 1-km Earth-crossers could have come from the Edgeworth-Kuiper belt into Jupiter- 
crossing orbits. Note that previous estimates of the number of bodies with diameter d>l km 
and 30<a<50 AU were larger: 10^° [12] and 10-*^^ [13] In the present paper we use the estimates 
from [2], but now include a much larger number of JCOs. Preliminary results were presented 
by Ipatov [14]-[16], who also discussed the formation of TNOs and asteroids, and in [17]. 

COLLISION PROBABILITIES IN THE MODEL OF FIXED ORBITAL ELE- 
MENTS 



Table 1. Characteristic collision times Tf (in Myr) of minor bodies with planets, coefficient 
k, and number Nf of simulated objects for the set of NEOs known in 2001. 
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In this section we estimate the probabilities of collisions of near-Earth objects with planets 
in the model of fixed orbital elements. As the actual collisions of migrating objects with terres- 
trial planets are rare, we use an approximation of random phases and orientations to estimate 
probabilities of collision for families of objects with similar orbital elements. We suppose that 
their semi-major axes a, eccentricities e and inclinations i are fixed, but the orientations of the 
orbits can vary. When a minor body collides with a planet at a distance R from the Sun, the 
characteristic time to collide, Tf, is a factor of A: = v/vc = yJlajR — 1 times that computed 
with an approximation of constant velocity, where v is the velocity at the point where the orbit 
of the body crosses the orbit of the planet, and Vc is the velocity for the same semi-major axis 
and a circular orbit. This coefficient k modifies the formulas obtained by Ipatov [1],[18] for char- 
acteristic collision and close encounter times of two objects moving around the Sun in crossing 
orbits. These formulas also depend on the synodic period and improve on Opik's formulas when 
the semi-major axes of the objects are close to each other. As an example, at e=0.7 and a=3.06 
AU, we have fc=2.26. 

Based on these formulas, we calculated probabilities {1/Tf) for ~ 1300 NEOs, including 343 
Venus-crossers, 756 Earth-crossers and 1197 Mars-crossers. The values of Tf (in Myr), k, and 
the number Nf of objects considered are presented in Table 1. We considered separately the 
Atens, Apollos, Amors, and several Jupiter-family comets (JFCs). The relatively small values 
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Table 2. Semi-major axes (in AU), eccentricities and inclinations of considered comets 
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of Tf for Atens and for all NEOs colliding with the Earth are due to several Atens with small 
inclinations discovered during the last three years. If we increase the inclination of the Aten 
object 2000 SG344 from i=0.1° to i=l°, then for colhsions with the Earth we find r=28 Myr 
and A;=0.84 for Atens and T=97 Myr and k=1.09 for NEOs. These times are much longer, and 
illustrate the importance of rare objects. 

INITIAL DATA 

As the migration of TNOs to Jupiter's orbit was investigated by several authors, we have 
made a series of simulations of the orbital evolution of JCOs under the gravitational influence 
of planets. We omitted the influence of Mercury (except for Comet 2P) and Pluto. The orbital 
evolution of about 9352 and 10301 JCOs with initial periods Pa<20 yr was integrated with 
the use of the Bulirsh-Stoer (BULSTO code [19]) and symplectic (RMVS3 code) methods, 
respectively. We used the integration package of Levison and Duncan [20] . 

In the first series of runs (denoted as nl) we calculated the evolution of 3100 JCOs moving in 
initial orbits close to those of 20 real comets with period 5<Pa<9 yr, and in the second series 
of runs (denoted as n2) we considered 7250 JCOs moving in initial orbits close to those of 10 
real comets (with numbers 77, 81, 82, 88, 90, 94, 96, 97, 110, 113) with period 5<Pa<15 yr. 
In other series of runs, initial orbits were close to those of a single comet (2P, 9P, lOP, 22P, 
28P, 39P, and 44P). In order to compare the orbital evolution of comets and asteroids, we also 
investigated the orbital evolution of asteroids initially moving in the 3:1 and 5:2 resonances with 
Jupiter. The number of objects in one run usually was <250. 

In most JCO cases the time r when perihelion was passed was varied with a step dr<l day 
(i.e., was varied with a step <0.2°). Near the r estimated from observations, we used smaller 
steps. In most JCO cases the range of initial values of r was less than several tens of days. 
For asteroids, we varied initial values of v and the longitude of the ascending node from to 
360°. The approximate values of initial orbital elements (a in AU, i in deg) are presented in 
Tabic 2. Wc initially integrated the orbits for r5>10 Myr. After 10 Myr we tested whether 
some of remaining objects could reach inside Jupiter's orbit; if so, the calculations were usually 
continued. Therefore the results for orbits crossing or inside Jupiter's orbit were the same as 
if the integrations had been carried to the entire lifetimes of the objects. For Comet 2P and 
resonant asteroids, we integrated until all objects were ejected into hyperbolic orbits or collided 
with the Sun. In some previous publications we have used smaller Tg, so these new data are 
more accurate. 

In our runs, planets were considered as material points so literal collisions did not occur. 
However, using the formulas of the previous section, and the orbital elements sampled with a 
500 yr step, we calculated the mean probability P of collisions. We define P as Pj:/N, where Ps 
is the probability for all N objects of a collision of an object with a planet during its lifetime, 
the mean time T=T^/N during which perihelion distance q of an object was less than the semi- 
major axis ttpi of the planet, the mean time T^, (in Kyr) spent in orbits with Q<4.2 AU, and the 
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Table 3. Mean probability P=W~^Pr of a collision of an object with a planet (Venus=V, Earth=E, 
Mars=M) during its lifetime, mean time T (in Kyr) during which q<api, Tc=T/P (in Gyr), mean time 
Tj (in Kyr) spent in Jupiter-crossing orbits, mean time (in Kyr) spent in orbits with Q<A.2 AU, 
and ratio r of times spent in Apollo and Amor orbits. Results from BULSTO code at 10^^<e<10^^ 
(marked as 10"^) and at e<ld~^'^ (marked as IQ-^^) and with RMVS3 code at integration step dg. 
The series of runs with a few excluded objects that had the largest probabilities of collision with the 
Earth are marked by *. 
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mean time Tj during which an object moved in Jupiter-crossing orbits. The values of Pr=lQ P, 
Tj, Td, and T are shown in Tables 3-4. Here r is the ratio of the total time interval when 
orbits are of Apollo type (a>l AU, q=a{l - e)<1.017 AU) at e<0.999 to that of Amor type 
(1.017<g<1.3 AU) and Tc=T/P (in Gyr). In almost all runs T was equal to the mean time in 
planet-crossing orbits and l/Tc was a probability of a collision per year (similar to 1/2/). 

In Table 5 for several objects that had large collision probabilities with the Earth, we present 
times (Myr) spent by these objects in orbits typical for inner-Earth objects (lEOs, (5<0.983 
AU), Aten (a<l AU and Q>0.983 AU), A12 (g<1.017 AU and l<a<2 AU), Apollo, and Amor 
objects, and also probabilities of collisions with Venus (pv), Earth (pe), and Mars (pm) during 
their lifetimes T^ (in Myr). Objects 44P and 113P presented in Table 5 were not included in the 
lines nl and n2 marked by * in Table 3, respectively, and objects 2P and lOP with BULSTO 
were not included in the lines for the corresponding series of runs in Table 3 with N=250 for 
2P and iV=2149 for lOP. 
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Table 4. Same as for Table 2, but for resonant asteroids. io=10° at eo=0.15, and io=5° at eo=0.05. 
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ORBITAL EVOLUTION OBTAINED BY DIRECT INTEGRATIONS 

Here and in Figs. 1, 2a-c, 3-5 we present the results obtained by the Bulirsh-Stoer method 
(BULSTO code [19]) with the integration step error less than eG[10-9-10-8], and in the next 
section we compare them with those of BULSTO with £<10-i2 and a symplectic method. 



Table 5. Times (Myr) spent by six objects in various orbits, and probabilities of collisions with 
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The results showed that most of the probability of collisions of former JCOs with the terrestrial 
planets is due to a small (^^0.1-1%) fraction that orbited for several Myr with aphelion Q<4.7 
AU. Some had typical asteroidal and NEC orbits and reached Q<3 AU for several Myr. Time 
variations in orbital elements of JCOs obtained by the BULSTO code are presented in Figs. 1, 
2a-b. Plots in Fig. 1 arc more typical than those in Fig. 2a-b, which were obtained for two 
JCOs with the highest probabilities of collisions with the terrestrial planets. Fig. 2c shows the 
plots for an asteroid from the 3:1 resonance with Jupiter. The results obtained by a symplectic 
code for two JCOs are presented in Fig. 2d-e. Large values of P for Mars in the nl runs were 
caused by a single object with a lifetime of 26 Myr. 

The total times for Earth-crossing objects were mainly due to a few tens of objects with high 
collision probabilities. Of the JCOs with initial orbits close to those of lOP and 2P, six and 
nine respectively moved into Apollo orbits with a<2 AU (A12 orbits) for at least 0.5 Myr each, 
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and five of them remained in sucli orbits for more than 5 Myr each. The contribution of ah the 
other objects to A12 orbits was smaller. Only one and two JCOs reached lEO and Aten orbits, 
respectively. 



Table 6. Times (in Myr) spent by N JCOs and asteroids during their lifetimes, with results for first 
50 Myr in [ ]. 
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One former JCO (Fig. 2a), which had an initial orbit close to that of lOP, moved in Aten 
orbits for 3.45 Myr, and the probability of its collision with the Earth from such orbits was 0.344 
(so rc=10 Myr was even smaller than the values of Tf presented in Table 1; i.e., this object had 
smaller e and i than typical observed Atcns), greater than that for the 9350 other simulated 
former JCOs during their lifetimes (0.17). It also moved for about 10 Myr in inner-Earth orbits 
before its collision with Venus, and during this time the probability i-V=0.655 of its collision 
with Venus was greater {Py^S for the time interval presented in Fig. 2a) than that for the 
9350 JCOs during their lifetimes (0.15). At t=0.12 Myr orbital elements of this object jumped 
considerably and the Tisserand parameter increased from J<3 to J>6, and J>10 during most 
of its lifetime. Another object (Fig. 2b) moved in highly eccentric Aten orbits for 83 Myr, 
and its lifetime before collision with the Sun was 352 Myr. Its probability of collisions with 
Earth, Venus and Mars during its lifetime was 0.172, 0.224, and 0.065, respectively. These two 
objects were not included in Table 3. Ipatov [21] obtained the migration of JCOs into lEO 
and Aten orbits using the approximate method of spheres of action for taking into account the 
gravitational interactions of bodies with planets. In the present paper we consider only the 
integration into the future. Ipatov and Hahn [9] integrated the evolution of Comet P/1996 R2 
both into the future and into the past, in this case the mean time Te during which a JCO was 
moving in Earth-crossing orbits is Tg = 5 x 10^ yr. The ratio Ps of the number of objects 
colliding with the Sun to the total number of escaped (collided or ejected) objects was less than 
0.015 for the considered runs (except for 2P). 

Ratio Ps of objects colliding with the Sun to those colliding with planets or ejected 
Series nl 9P lOP 22P 28P 39P 44P 
Ps 0.0005 0.014 0.002 0.007 0.004 

Some former JCOs spent a long time in the 3:1 resonance with Jupiter (Fig. la-b) and with 
2<a<2.6 AU. Other objects reached Mars-crossing orbits for long times. We conclude that 

JCOs can supply bodies to the regions which are considered by many scientists [3] to belong to 
the main sources of NEOs, and that those rare objects that make transitions to typical NEO 
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Fig. 1. Time variations in a, e, q, Q, sin(i) for a former JCO in initial orbit close to that of Comet 
lOP (a-f), or Comet 2P (g-h). Results from BULSTO code with s ~ IQ-^ - IQ-^. 
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Fig. 2. Time variations in a, e, q, Q, and i for a former JCO in initial orbit close to that of Comet 
lOP (a), 2P (b, e), 9P (d), or an asteroid at the 3/1 resonance with Jupiter (c). For (a) at i<0.123 
Myr Q>a>1.5 AU. Results from BULSTO code with e ~ 10~^ - 10~^ (a-c) and by a symplectic 
method with ^^=30 days (d) and with dg—lO days (e). 
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orbits dominate the statistics. Only computations with very large numbers of objects can hope 
to reach accurate conclusions on collision probabilities with the terrestrial planets. 

In Fig. 3 we present the time in Myr during which objects had semi-major axes in an interval 
with a width of 0.005 AU (Figs. 3a-b) or 0.1 AU (Figs. 3c). At 3.3 AU (the 2:1 resonance with 
Jupiter) there is a gap for asteroids that migrated from the 5:2 resonance and for former JCOs 
(except 2P). 

For the nl data set, r7=0.12 Myr and, while moving in Jupiter-crossing orbits, objects had 
orbital periods Pa<10, 10<P„<20, 20<Pa<50, 50<Pa<200 yr for 11%, 21%, 21%, and 17% of 
Tj, respectively. Therefore, there are three times as many JCOs as Jupiter- family comets (for 
which Pa<20 yr). We also found that some JCOs, after residing in orbits with aphelia deep 
inside Jupiter's orbit, transfer for tens of Myr to the trans-Neptunian region, either in low or 
high eccentricity orbits. We conclude that some of the main belt asteroids may reach typical 
TNO orbits, and then become scattered-disk objects having high eccentricities, and vice versa. 
The fraction of objects from the 5:2 resonance that collided with the Earth was only 1/6 of that 
for the 3:1 resonance. Only a small fraction of the asteroids from the 5:2 resonance reached a<2 
AU (Fig. 3b). 

The distributions of migrating former JCOs (2P and lOP) and resonant asteroids in a and 
e (left) and in a and i (right) are presented in Fig. 4-5. For each picture we considered 250 
migrating objects (288 for Fig. 5), 100 intervals for a, and about the same number of intervals for 
e and i. Different designations correspond to different numbers n of orbital elements (calculated 
with a step of 500 yr) in one bin (in Fig. 4 '<=' means <). Similar plots for 39P runs were 
presented in [17]. All the former JCOs reached low eccentricity orbits very rarely with 2<a<3.5 
AU and ll<a<28 AU. There were many positions of objects when their perihelia were close 
to a semi-major axis of a giant planet, mainly of Jupiter (Fig. 4a). Note that Ozernoy et al. 
[22] considered the migration of Neptune-crossers and found that the main concentrations of 
perihelia were near Neptune's orbit. 

W. Bottke pointed out that H. Levison showed that it is difficult to detect solar collisions 
in any numerical integrator, so he removed objects with g<gmin- The results presented above 
were obtained considering collisions with the Sun, but we also investigated what happens if we 
consider qrahi equal to ks radii rs of the Sun. For ks=2, some results are presented in Fig. 4d-e, 
5b. The only difference with the runs that considered collisions with the Sun is that for those 
runs for series 2P and lOP and for the 3:1 resonance, some objects reached 90°<i<180° (mainly 
with 2<a<3.5 AU) (Fig. 4b-c, 5a). For ks=2 there were no comets with i>90° and there were 
only a few orbits of asteroids with i>90° (Fig. 4d-e, 5b). The consideration of qinin at ks=3 did 
not influence the collision probabilities with the terrestrial planets or getting orbits with a<2 
AU. For example, with BULSTO for the two objects with the largest collision probabilities, the 
time spent in orbits with a<2 AU decreased by only 0.3% for 2P at ks=S and was the same for 
lOP at ks=lO. 

COMPARISON OF ORBIT INTEGRATORS 

To determine the effect of the choice of orbit integrators and convergence criteria, we made 
additional runs with BULSTO at e=10~^^ and £=10~^^ and with a symplectic integrator. The 
orbital evolution of 9551 JCOs was computed with the RMVS3 code. For the symplectic method 
we usually used an integration step dg of 3, 10, and 30 days. For series n2 and 44P we took 
different values of dg between 5 and 10 days. 

We find that, exclusive for the case of close encounters with the Sun, the differences be- 
tween integrator choices (with ds<lO days) are comparable to the differences between runs with 
slightly different initial conditions. Our interpretation is that 1) very small numbers of particles 
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Fig. 3. Distribution of 7852 migrating JCOs (a, c) and 288 resonant asteroids at eo=0.15 and 
io=10° (b) with their semi-major axes. The curves plotted in (c) at a=40 AU are (top-to-bottom) 
for sum, lOP, nl, 39P, 22P, 9P, 28P, and 2P (series n2 and 44P are not included in the figure). For 
Figs, (a) and (c), designations are the same. Results from BULSTO code with e ~ 10~^ — 10~^. 
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Fig. 4. Distribution of migrating objects in semi-nnajor axes, eccentricities, and inclinations for objects 
in initial orbits close to that of lOP (a-b,d), 2P (c,e), and BULSTO code with e ~ 10"^ - 10"^. 
For (d-e) it was considered that an object disappeared when perihelion distance became less than 2 
radii of the Sun. In other cases (a-c) objects disappered when they collided with the Sun. 
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Fig. 5. Distribution of migrating objects in semi-major axes and inclinations for objects in initial orbits 
at the 3:1 resonance with Jupiter with eo=0.15 and io—10°. BULSTO code with e ~ 10~^ — 10"^. 
For (a) objects disappered when they collided with the Sun. For (b) it was considered that an object 
disappeared when perihelion distance became less than 2 radii of the Sun. 



contribute most of the collision probabilities with the terrestrial planets, 2) runs with larger 
numbers of particles are more reliable, and 3) small differences in initial conditions or in the 
errors of the orbit integrators modify the trajectories substantially, especially for those particles 
making major changes in their orbits due to close encounters or resonances. We eonelude that 
for the purposes of this paper, the various choices of orbit integrator are sufficiently equivalent. 

To illustrate these points. Tables 3-4 present the results obtained by BULSTO with e<10~^^ 
and the symplectic method with ds<10 days. Most of the results obtained with these values 
of £ and ds are statistically similar to those obtained for 10~®<e<10~^. For example, a few 
objects spent millions of years in Earth-crossing orbits inside Jupiter's orbit (Figs. 1-2), and 
their probabilities of collisions with the Earth were thousands of times greater than for more 
typical objects. For series nl with ds<lO^, the probability of a collision with Earth for one 
object with initial orbit close to that of Comet 44P was 88.3% of the total probability for 1200 
objects from this series, and the total probability for 1198 objects was only 4%. This object and 
the object presented in Fig. 2e and in the third line of Table 5 were not included in Table 3 with 
iV=1199 for nl and with Ar=250 for 2P, respectively. For the 3:1 resonance with ^5=10 days, 
142 objects spent 140 and 84.5 Myr in lEO and Aten orbits, respectively, even longer than for 
e~10~^-10^^. Additionally, up to 40 Myr and 20 Myr were spent in such orbits by two other 
objects which had estimated probabilities of collisions with the terrestrial planets greater than 
1. For the 2P runs with £:<10~^^ and A/"=100, the calculated objects spent 5.4 Myr in Apollo 
orbits with a<2 AU. 

The values of Pr are usually of the same order of magnitude for different methods (see Tables 
3-4), and the difference between the data is comparable to the differences between different runs 
belonging to a series. For Earth and Venus, the values of Pr are about 1-4 for Comets 9P, 
22P, 28P and 39P; for Comet 44P they are about 4-6. For 28P and 39P with the symplectic 
method Pr is about twice that for BULSTO. For lOP the values of Pr are about 20-40 and are 
several times larger than for the above series, and for 2P they exceed 100 and are several times 
larger than for lOP. With the nl and n2 runs, for Earth Pj.>4 and Pr.>10, respectively. The 
ratio of Pr to the mass of the planet was typically several times larger for Mars than for Earth 
and Venus. The main difference in Pr was found for the 3:1 resonance. In this case greater 
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values of P,,. were obtained for (is=10 days than for BULSTO. As noted above, a few exceptional 
objects dominated the probabilities, and for the 3:1 resonance two objects, which had collision 
probabilities greater than unity for the terrestrial planets, were not included in Table 4. These 
two objects can increase the total value of Pr for Earth by a factor of several. 

The mean time T^^ (in Kyr) spent in orbits with Q<4.2 AU can differ by three orders of 
magnitude for different series of runs (Tables 3-4) . For most runs (except for 2P and asteroids) 
the number of objects which got Q<4:.7 AU was several times larger than that for Q<4:.2 AU. 

For symplectic runs with ^^=30 days for most of the objects we got results similar to those 
with ds<10 days, but about 0.1% of the objects reached Earth-crossing orbits with a<2 AU for 
several tens of Myr (e.g., Fig. 2d) and even lEO orbits. These few bodies increased the mean 
value of P by a factor of more than 10. With ds=30 days, four objects from the runs nl, 9P, 
lOP had a probability of collisions with the terrestrial planets greater than 1 for each, and for 
2P there were 21 such objects among 251 considered. For resonant asteroids, we also obtained 
much larger values than those for BULSTO for P and T for RMVS3 with ds=30 days, and 
similarly for the 3:1 resonance even with ds=W days. For this resonance it may be better to 
use ds<10 days. Probably, the results of symplectic runs with cZs=30 days can be considered as 
such migration that includes some nongravitational forces. 

In the case of close encounters with the Sun (Comet 2P and resonant asteroids), the probability 
Ps of collisions with the Sun during lifetimes of objects was larger for RMVS3 than for BULSTO, 
and for 10~^^<e<10~^^ it was greater than for 10~^<e<10^® (P5'=0.75 for the 3:1 resonance 
with ds=3 days). This probability is presented in Table 7 for several runs. 



Table 7. Probability of collisions with the Sun (for asteroids eo=015 and io—lO°). 
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For Comet 2P the values of Tj were much smaller for RMVS3 than those for BULSTO and 
they were smaller for smaller e; for other runs these values do not depend much on the method. 
In our opinion, the most reliable values of Tj were obtained with 10~^^<£<10~^^. In the direct 
integrations reported by Valsecchi et al. [23], 13 of the 21 objects fell into the Sun, so their 
value of ^5=0. 62 is in accordance with our results obtained by BULSTO; it is less than that 
for e=10^^^, but greater than for e=10~^. Note that even for different Ps the data presented 
in Tables 2-3 usually are similar. As we did not calculate collision probabilities of objects with 
planets by direct integrations, but instead calculated them with the random phase approximation 
from the orbital elements, we need not make integrations with extremely high accuracy. We 
showed [24] that for BULSTO the integrals of motion were conserved better and the plots of 
orbital elements for closely separated values of e were closer to one another with 10~^<e<10~^. 
The smaller the value of e, the more integrations steps are required, so £:<10~^^ for large time 
intervals are not necessarily better than those for 10^^<e<10^*^. Small e is clearly necessary 
for close encounters. Therefore we made most of our BULSTO runs with 10~^<e<10^^. We 
found [1],[9] that former JCOs reached resonances more often for BULSTO than for RMVS3 
with ds=SO days. For a symplectic method it is better to use smaller dg at a smaller distance 



14 



R from the Sun, but in some runs R can vary considerably during the evolution. The choice of 
dg depends on the smallest values of R, so a symplectic method with a constant ds may not be 
effective when R is very different for different objects considered. 

MIGRATION FROM BEYOND JUPITER TO THE EARTH 

The fraction Ptnj of TNOs reaching Jupiter's orbit under the influence of the giant planets 
in 1 Gyr is 0.8-1.7% [6]. As the mutual gravitational influence of TNOs can play a larger role in 
variations of their orbital elements than collisions [2], we considered the upper value of Ptnj- 
Using the total of 5 x 10^ 1-km TNOs with 30<a<50 AU, and assuming that the mean time 
for a body to move in a Jupiter-crossing orbit is 0.12 Myr, we find that about Ajo= 10^ 1-km 
former TNOs are now Jupiter-crossers, and 3000 are Jupiter-family comets. With the total 
times spent by Nj simulated JCOs in Apollo orbits we can estimate the number of 1-km former 
TNOs now moving in such orbits using the following formula: N Apollo = ^Jo ' t Apollo/ J ■ ij), 
where t Apollo is the total time during which Nj former JCOs moved in Apollo orbits, and Nj ■ tj 
is the total time during which Nj JCOs moved in Jupiter-crossing orbits. Similar formulas can 
be considered for other types of orbits. Based on nl and n2 series of runs, we obtain that there 
are about 700-900 Amors and 2000-3000 Apollos (the last numbers include very eccentric orbits) 
which came from the Edgeworth-Kuiper belt and have diameters greater than 1 km. Even if the 
number of Apollo objects is an order of magnitude smaller than the above value, it may still be 
comparable to the real number (750) of 1-km Earth-crossing objects (half of them are in orbits 
with a<2 AU), although the latter number does not include those in highly eccentric orbits. 

The ratio k2 of the number of Apollos with a<2 AU to the number of all Apollos was very 
different for different series of runs. It was 0.4 for nl runs, but the total time spent by 3100 
considered objects in orbits with a<2 AU was due mainly only to one object. For n2 series with 
RMVS3, A;2=0.2, but practically all time spent by 6250 objects in orbits with a<2 AU was also 
due to one object. For all JCOs considered, k2 was about 0.5, but the total time spent in orbits 
with a<2 AU was due mainly to the objects from 2P series of runs. 

The values of the characteristic time (usually for the collision of a former JCO or a resonant 
asteroid with a planet (see Tables 3-4) are greater than the values of Tf for NEOs in Table 1, 
so we expect that the mean inclinations and eccentricities of unobserved NEOs are greater than 
those for the NEOs that are already known. Similar results were found in [25]. On average, 
the values of Tc for our nl and n2 series and for most of our simulated JCOs were not greater 
than those for our calculated asteroids, and migrating Earth-crossing objects had similar e and 
i for both former JCOs and resonant asteroids. Former JCOs, which move in Earth-crossing 
orbits for more than 1 Myr, while moving in such orbits, usually had larger P and smaller e 
and i (sometimes similar to those of the observed NEOs, see Figs. 1-2). It is easier to observe 
orbits with smaller values of e and i, and probably, many of the NEOs moving in orbits with 
large values of e and i have not yet been discovered. About 1% of the observed Apollos cross 
Jupiter's orbit, and an additional 1% of Apollos have aphelia between 4.7-4.8 AU, but these 
Jupiter-crossers are far from the Earth most of time, so their actual fraction of ECOs is greater 
than for observed ECOs. The fraction of Earth-crossers among observed Jupiter-family comets 
is about 10%. This is a little more than T/Tj for our nl runs, but less than for n2 runs. For 
our former resonant asteroids, Tj is relatively large (^0.2 Myr), and such asteroids can reach 
cometary orbits. 

Comets are estimated to be active for Tact'~^10^-10'^ yr. Tact is smaller for closer encounters 
with the Sun [5], so for Comet 2P it is smaller than for other JFCs. Some former comets can 
move for tens or even hundreds of Myr in NEO orbits, so the number of extinct comets can 
exceed the number of active comets by several orders of magnitude. The mean time spent by 
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Encke-type objects in Earth-crossing orbits is >0.4 Myr (even for q„im)- This time corresponds 
to >40-400 extinct comets of this type. Note that the diameter of Comet 2P is about 5-10 km, 
so the number of smaller extinct comets can be much larger. 

The above estimates of the number of NEOs are approximate. For example, it is possible that 
the number of 1-km TNOs is several times smaller than 5 x 10^, while some scientists estimated 
that this number can be up to 10^"*^ [13]. Also, the fraction of TNOs that have migrated towards 
the Earth might be smaller. On the other hand, the above number of TNOs was estimated for 
a < 50 AU, and TNOs from more distant regions can also migrate inward. Probably, the Oort 
cloud could also supply .jTipitcr-faniily comets. According to Asher et al. [26], the rate of a 
cometary object decoupling from the Jupiter vicinity and transferring to an NEO-like orbit is 
increased by a factor of 4 or 5 due to nongravitational effects (see also [27]). This would result 
in larger values of Pr and T than those shown in Table 3. 

Our estimates show that, in principle, the trans-Neptunian belt can provide a significant 
portion of the Earth-crossing objects, although many NEOs clearly came from the main asteroid 
belt. Many former Jupiter-family comets can have orbits typical of asteroids, and collide with 
the Earth from typical NEO orbits. It may be possible to explore former TNOs near the Earth's 
orbit without sending spacecraft to the trans-Neptunian region. 

Based on the estimated collision probability P = 4 x 10~^ wc find that 1-km former TNOs 
now collide with the Earth once in 3 Myr. This value of P is smaller than that for our nl, and 
especially than for n2, lOP and 2P runs. Assuming the total mass of planetesimals that ever 
crossed Jupiter's orbit is ~lOOm0, where is the mass of the Earth [1],[28], we conclude that 
the total mass of bodies that impacted the Earth is 4 x lO"^^®. If ices comprised only half of 
this mass, then the total mass of ices Mice that were delivered to the Earth from the feeding 
zone of the giant planets is about the mass of the terrestrial oceans (~ 2 x lO^^m®). 

The calculated probabilities of collisions of objects with planets show that the fraction of 
the mass of the planet delivered by short-period comets can be greater for Mars and Venus 
than for the Earth (Table 3). This larger mass fraction would result in relatively large ancient 
oceans on Mars and Venus. On the other hand, there is the deuterium/hydrogen paradox of 
Earth's oceans, as the D/H ratio is different for oceans and comets. Pavlov et al. [29] suggested 
that solar wind-implanted hydrogen on interplanetary dust particles could provide the necessary 
low-D/H component of Earth's water inventory, and Delsemme [30] considered that most of the 
seawater was brought by the comets that originated in Jupiter's zone, where steam from the 
inner solar system condensed onto icy interstellar grains before they accreted into larger bodies. 

Our estimate of the migration of water to the early Earth is in accordance with [31], but 
is greater than those of Morbidelli et al. [32] and Levison et al. [33]. The latter obtained 
smaller values of Mice, and we suspect that this is because they did not take into account 
the migration of bodies into orbits with Q<4.5 AU. Perhaps this was because they modeled 
a relatively small number of objects, and Levison et al. [33] did not take into account the 
influence of the terrestrial planets. In our runs the probability of a collision of a single object 
with a terrestrial planet could be much greater than the total probability of thousands of other 
objects, so the statistics are dominated by rare occurrences that might not appear in smaller 
simulations. The mean probabilities of collisions can differ by orders of magnitude for different 
JCOs. Other scientists considered other initial objects and smaller numbers of Jupiter-crossing 
objects, and did not find decoupling from Jupiter, which is a rare event. We believe there is 
no contradiction between our present results and the smaller migration of former JCOs to the 
near-Earth space that was obtained in earlier work, including our own papers (e.g. [9]), where 
we used the same integration package. 
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Prom measured albedos, Fernandez et al. [34] concluded that the fraction of extinct comets 
among NEOs and unusual asteroids is significant (at least 9% are candidates). The idea that 
there may be many extinct comets among NEOs was considered by several scientists. Rickman 
et al. [35] believed that comets played an important and perhaps even dominant role among all 
km-size Earth impactors. In their opinion, dark spectral classes that might include the ex-comets 
are severely underrepresented (see also [10]). Our runs showed that if one observes former comets 
in NEO orbits, then it is probable that they have already moved in such orbits for millions (or 
at least hundreds of thousands) years, and only a few of them have been in such orbits for short 
times (a few thousand years). Some former comets that have moved in typical NEO orbits for 
millions or even hundreds of millions of years, and might have had multiple close encounters with 
the Sun (some of these encounters can be very close to the Sun, e.g. in the case of Comet 2P at 
t>0.05 Myr), could have lost their mantles, which causes their low albedo, and so change their 
albedo (for most observed NEOs, the albedo is greater than that for comets [34]) and would 
look like typical asteroids or some of them could disintegrate into mini-comets. Typical comets 
have larger rotation periods than typical NEOs [36]-[37], but, while losing considerable portions 
of their masses, extinct comets can decrease these periods. For better estimates of the portion 
of extinct comets among NEOs we will need orbit integrations for many more TNOs and JCOs, 
and wider analysis of observations and craters. 

CONCLUSIONS 

Some Jupiter-family comets can reach typical NEO orbits and remain there for millions of 
years. While the probability of such events is small (about ~0.1%), nevertheless the majority 
of collisions of former JCOs with the terrestrial planets are due to such objects. Most former 
TNOs that have typical NEO orbits moved in such orbits for millions of years (if they did not 
disintegrate into mini-comets), so during most of this time there were extinct comets. From the 
dynamical point of view there could be many extinct comets among the NEOs. The amount 
of water delivered to the Earth during planet formation could be about the mass of the Earth 
oceans. 
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